Dynamical spin structure factors of quantum spin nematic states 
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Dynamical spin structure factors of quantum spin nematic phase are studied using the large- 
N loop expansion in the spin-i square-lattice J1-J2 model with ferromagnetic Ji and competing 
antiferromagnetic J2. As the starting mean- field state, we employ a spin-triplet pairing state of 
spinon fields, called Z2 planar state, which was found as a large-A'' saddle point solution of the 
fermionic mean-field analysis [R. Shindou and T. Momoi, Phys. Rev. B 80, 064410 (2009)]. Using 
the standard large- A'^ expansion, we take into account the fluctuation within 1-loop level, which is 
equivalent to the random phase approximation. The spin structure factors thus obtained signifies the 
existence of the gapless fc-linear spin- wave modes at q = (0,0) and that at q = (0,7r). The spectral 
weight of the former gapless modes vanish as linear functions of the momentum, while the latter 
mode has essentially no spectral weight in the dynamical spin structure factors. We also observed 
that the first gapped mode at q = (7r,7r), which can be observed in the longitudinal dynamical 
spin-structure factor, is a 'Higgs boson' associated with the Z2 state. Namely, with decreasing J2, 
the mass of the 'Higgs boson' vanishes at a critical value J2/J1 ~ Jc,2, so that the state reduces 
to the (7(1) planar state. We found that, when this happens, other gapped magnetic modes at 
q = (7r,7r) simultaneously exhibit instabilities. As a result, the U{1) planar phase not only breaks 
the translational symmetries of the square lattice, but also it breaks a staggered 17(1) spin-rotational 
symmetry which is possessed by the Z2 planar state. 

PACS numbers: 



I. INTRODUCTION 

Frustrated magnets are Mott insulators where com- 
peting interactions among localized quantum spins bring 
about a large degeneracy in the ground state energet- 
ics. In a certain circumstance, such a frustrated spin 
system lifts this degeneracy quantum-mechanically, only 
to choose as its ground state a liquid-like magnetic state 
of matter, dubbed as a quantum spin liquid ^i, where the 
system does not exhibit any ordering down to the zero 
temperature. 

The representative of quantum spin liquids is resonat- 
ing valence bond (RVB) state in S* = i quantum spin sys- 
tems;^ whose building-block is the singlet valence bond, 
a pair of two S — ^ spins forming the spin-singlet state. 
The wavefunction is given by an equally-weighted super- 
position of different spin-singlet dimer covering configu- 
rations, so that, unlike valence bond solid state or Neel 
state, the state preserves the lattice-translational symme- 
try. Owing to this fluid-like character of the ground-state 
wavefunction, quantum spin liquids have various exotic 
low-energy excitations, including emergent gauge fields, 
fractionalized magnetic excitations called spinous, and 
the topological degeneracy associated with this fraction- 
alizatioufir— 

A new state of matter recently discussed in the lo- 
calized spin systems is a quantum spin nematic state, 
which can be regarded as a quantum-spin analogue of 
nematic liquid crystals. It neither possesses spin order, 
i.e. sublattice magnetization, nor crystalline solid-like 
structure in spin degrees of freedom, but, unlike spin- 



rotational symmetric quantum spin liquids, it exhibits 
the spin nematic order measured by the rank-2 traceless 
tensor spin operatori^ Kjrn,tii' = Sj^f,Sm,i' -^{Sj ■ Sm) 
for /i, = 1,2,3. The tensor is comprised by two dis- 
tinct spin operators defined on different sites, usually 
neighboring two sites, so that the tensor is introduced 
on the bond. The state is thereby called as the 'bond- 
type' spin nematic state. Depending on the symme- 
try under the exchange of spin indices, fi ^ v, quan- 
tum spin nematics can be classified into two categories)^ 
the chiral type (p-nematic) and the non-chiral type (n- 
nematic). The order parameter of the former one is 
the vector chirality defined as {Sj x Sm), while that of 
the latter is the symmetric part of the quadratic tensor. 



{Qj 



A p-nematic state was first introduced as a helical 
quantum spin liquid in a frustrated antiferromagnetic 
Heisenberg model by Chandra, Coleman and Larkin.— 
The state essentially acquires a linear combination of the 
spin- triplet valence bond and singlet valence bond on the 
same bond, i.e. | fD+e'"^! it): which induces a finite vec- 
tor chirality on the bond. Meanwhile, possible realization 
of n-nematic states has been investigated recently in var- 
ious quantum frustrated ferromagnetic Heisenberg mod- 
els,—"— where dynamics of (purely) spin- triplet valence 
bonds brings about the long range ordering of the sym- 
metric order parameter {Qjm,iiv)- Both of these bond- 
type spin nematic states preserve the translational sym- 
metries of the original lattices, so that they can be re- 
garded as a 'cousin' of symmetric quantum spin liquids, 
sharing many of their exotic aspects. 

Motivated by these recent findings, the present authors 
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FIG. 1: Square lattice J1-J2 frustrated ferromagnetic model, 
where the nearest neighbor exchange Ji is ferromagnetic and 
the next nearest neighbor exchange J2 is antiferromagnetic. 



n-flux 


Zj planar 


U(l) planar 


Ferro- 


states 
■* 1 


state 
1 


state 
-i 1 


state 



■J 2 / 1 ^C,l ^C,2 ^c,3 

FIG. 2: Phase diagram of the square lattice J1-J2 model in 
the large- A*' limit. The nearest neighbor exchange Ji is ferro- 
magnetic and the next nearest neighbor exchange J2 is antifer- 
romagnetic. Energetics of the mean-field solutions conclude 
that Jc,i = 1.325, Jc,2 = 1.0448 and J^.s = 1.02. 



formulated the mean-field theory of bond-type spin ne- 
matic states in spin-i spin systems in terms of fermion 
representation, where the spin-triplet pairing of spinon 
fields describes the spin-triplet valence bond,^*^ It turns 
out that the so-called d- vector associated with this triplet 
pairing function mimics the director vector in the n- 
nematics, while the product between the d-vector and 
singlet pairing amplitude defines the vector chirality as- 
sociated with the p-nematics. Using this formulation, we 
derived the mean-field phase diagram of the square lattice 
J1-J2 model with ferromagnetic (FM) Ji and competing 
antiferromagnetic (AF) J2, thereby finding that a cer- 
tain spin triplet pairing state, dubbed as the Z2 planar 
state, is most energetically favorable when competition 
between the FM Ji and the AF J2 interactions is strong 
(see Fig. EJ. The variational Monte Carlo (VMC) study 
based on this mean-field analysis further indicated that, 
when projected onto the physical spin Hilbert space, the 
BCS wavefunction constructed from the Z2 planar state 
achieves the best optimal energy, compared with other 
competing states such as the ferromagnetic state and 
the coUinear antiferromagnetic state. ^'^ This projected Z2 
planar state possesses the same 'd-wave' spatial configu- 
ration of the bond quadrupole order— as suggested by 
the previous exact diagonalization studyi^ 

In this paper, we theoretically investigate the nature 
of low-energy magnetic excitations in the quantum spin 
nematic state. We calculate the dynamical magnetic 
properties of the Z2 planar state, so as to give relevant 
physical characterizations to the bond-type spin nematic 
phase. Specifically, using a standard large-iV loop ex- 
pansion, we take into account the fluctuation around the 
mean-field state within the one-loop level. This treat- 



ment essentially corresponds to the random phase ap- 
proximation (RPA). The spin structure factor thus ob- 
tained has two aspects; spin-liquid like character and 
spontaneously symmetry breaking (SSB) phase charac- 
ter. The former feature manifests itself as the Stoner 
continuum of the individual excitations of gapped free 
spinous. On the one hand, the SSB phase character is 
represented by the gapless spin-wave modes, which have 
linear energy dispersions and whose spectral weight van- 
ishes as a linear function of the momentum near the gap- 
less point. 

When the antiferromagnetic coupling J2 decreases, the 
mean-field solution transforms from the Z2 planar state 
to the U{1) planar state at J2/J1 = Jc,2 (see Fig. 

Namely, previous field-theoretical analysis suggested 
that, some of the gapped gauge bosons in the Z2 pla- 
nar state become gapless at the transition point, only to 
constitute a compact QED (quantum electrodynamics) 
action for J2/J1 < Jc,2, where the space-time instanton 
effect associated with this effective action introduces a 
strong confining potential between 'free' gapped spinous. 
Throughout the RPA calculation in this paper, we found 
that the first gapped modes at q = {t^jT^) are actually 
these gapped gauge bosons, whose finite mass therefore 
quantifies the stability of the Z2 planar state against the 
confinement effect. We also found that, when the Z2 
state is transformed into the U{1) state, two other mag- 
netic modes at the same zone boundary point simultane- 
ously exhibit the instability. As a result, the U{1) planar 
phase breaks the translational symmetries and also a cer- 
tain type of the spm-rotational symmetry. 

The organization of this paper is as follows: In the 
next section, we first introduce a 'general-A^' quantum 
frustrated ferromagnetic model, whose large- iV limit pos- 
sesses our previous mean-field solutions as the 'exact' 
ground states and whose N — 1 case safely reproduces 
the usual quantum spin-i model. In Sec. Ill, we will 
describe the -^-expansion calculation for the dynamical 
spin-spin correlation functions. In Sec. IV, we show the 
calculated dynamical spin structure factors, both longi- 
tudinal and transverse, and discuss their characteristic 
features and its physical implications. We also discuss 
about the nature of the U{1) planar state here. Section V 
is devoted to the summary and discussion. 



II. QUANTUM FRUSTRATED 
FERROMAGNETIC MODEL 



In this section, we will briefly review the Z2 planar 
state as a saddle point solution in the large- limit of 
a 'general-A^' quantum frustrated ferromagnetic model. 
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The large- iV model is given hyi^ 



n 



J ^ 

{j .,m) a, 6—1 



N 



J2 E-^f •^-+E'^r-'^r' w 

{{3,m)) a,b j ,a 

where {j,m) {{{j,m))) runs over all nearest- neighbor 
(2nd-neighbor) bonds on the square lattice and 
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(2) 



Here f'^'^^ is a fermion creation operator with spin 
a =t or \. and flavor a — \, . . . , N . In this paper, we con- 
sider the case Ji is ferromagnetic and J2 is antiferromag- 
netic, i.e. Ji, J2 > 0. We have introduced an external 
magnetic field h"-"" to calculate the spin-spin correlation 
function. The physical spin Hilbert space satisfies the 
constraint 



AT 



a=l 



N 

a=l 



on every site j. When iV = 1, Eq. ^ in the phys- 
ical Hilbert space reduces to the usual J1-J2 quantum 
Heisenberg model with spin-i. We regard N to be large 
when we perform the 1/iV loop-expansion. 

An equivalent statistical-mechanics problem at tem- 
perature j3~^ can be formulated in terms of the path- 
integral, 

Z[h] = 



exp 




■T S (-1 



dTC[h,U'"\U'",ar 



2+tr[*f[/j;;^,,*;;^aJ] 



(4) 



where the fermion field is written in the 2x2 matrix form 

(5) 

and the auxiliary fields are written as 



f "t fa 
faf _ ra 

Jj^i Jit 



^ jm 



X 



jm 'ijm 



u: 



tri 



FT"* r)* 

jm,fi jm.fj. 



(6) 



(/i = 1,2,3). The trace denoted by the symbol "tr" is 
taken over 2x2 matrices such as ^E"^ and the Pauli ma- 
trices cr^. Integrating over the temporal gauge fields aj -r 
strictly imposes the local constraints given by Eq. ([3]) 
on every site and time. In this formulation, the ferro- 
magnetic exchange interaction was decoupled in terms of 
spin-triplet pairing fields and spin-triplet hopping fields 

while antiferromagnetic one was decomposed into the 
spin-singlet pairing and hopping fields 

Vjl^-^ifjA'^^Uliflj), Xjl= {fljl,a)- (8) 

A Gaussian-integral over these auxiliary fields reproduces 
the original Hamiltonian given in Eq. ([TJ. 

We can formally rewrite the effective action as 



drC = NSi 

(9) 



fc,n,a 



introducing the Nambu representation of the fermion 
field 

fk,n,a = ( /fc,n,t fk,n,l f-k,-n,t f-k,-n,i ) i^^) 



with 



rat 



1 ri^ 

= Y dTe''=-^'+'""V"U^), (11) 



where N\ denotes the total number of lattice site and 
U}n = {2n + l)TTp-\ In Eq. ©, 



{j.m.) 

+ X E {\Xjm\ + \Vjrr.\)} (12) 

and denotes the inverse of a single-particle Green 
function of the fermion field. Integral over the ^ (/) 



4 



fields in Eq. Q leads to the following partition function, 



Z[h] 



J VW'"'VU'"Var exp [ - NS[h, 



tri 



S = Si + Su 

a=l 



(13) 
(14) 

(15) 



where the Green function G is diagonal in the flavor in- 
dex. The trace of lnG~^ is taken over the momentum 
(fc), the Matsubara frequency (w„), spin index (cr), and 
particle-hole index. 

iSii as well as Si are at most on the order of unit in 
large- iV limit, and the same is for S. Hence, in this limit, 
we can evaluate the partition function, substituting the 
saddle point solution of the action given by 



6S 



6S 



SS 



0. 



(16) 



The present authors previously investigated various local 
minima of the action at zero magnetic field, /i = 0, by 
assuming that C/'" and Ot- thus determined are tem- 
porally uniform and also preserve the translational sym- 
metries of the bsquare lattice. Throughout this study, 
we found that, in the intermediate coupling region, the 
saddle point solution acquires a finite spin-triplet pairing 
on every ferromagnetic bonds, while singlet pairings on 
the antiferromagnetic bonds. 

In particular, when the next nearest neighbor antifer- 
romagnetic exchange interaction supports the '7r-fiux'- 
type singlet pairings, a coplanar configuration of the spin- 
triplet d- vectors on the ferromagnetic interactions real- 
izes the 'best' mean-field energy among the others This 
mean-field solution, which we call 'Z2 planar state', is 
given by 



(17) 



with j = Uxjy), = (1,0) and By = (0,1). The in- 
variant gauge group of this mean-field state dictates that 
all the gauge excitations around this mean-field solution 
have finite gap [are not required to be gapless by the 
local SU{2) gauge symmetry]. The state has the same 
spin-triplet pairing function as a 'planar' type superfluid 
B-phase '^He: it is thereby dubbed as the 'Z2 planar 
state'. Though the coplanar ordering of the d- vectors 
breaks the SU (2) spin-rotational symmetries completely, 
this ordering also endows the state with a staggered U{1) 
spin-rotational symmetry. That is, the state defined by 
Eq. (|17p is invariant under the following rotation in the 
spin-space. 



, (18) 



for any 9, where the rotational axis is perpendicular to 
the coplanar plane. 



The variational Monte Carlo study further indicates 
that the projected BCS wavef unction constructed from 
this Z2 planar state achieves the best optimal energy in 
the range 0.42 Ji < J2 < 0.57Ji in intermediate coupling 
regime, when compared with energies of other compet- 
ing phases such as ferromagnetic state and coUinear Neel 
stateji^ Moreover, the projected BCS wavefunction be- 
longs to the same space group (including its irreducible 
representation)^ as that of a bond-type spin nematic 
phase suggested by the exact diagonalization (ED) study 
in the same parameter region. The spin-spin correlation 
function calculated with respect to this projected BCS 
wavefunction exhibits a similar behavior— as those ob- 
tained from the ED studies up to 40 sitesJ^ Observing 
these consistencies with the previous ED results includ- 
ing its energetics, we regard that this projected Z2 pla- 
nar phase is indeed realized as a bond-type spin nematic 
phase within a certain range of the intermediate coupling 
regime of the present J1-J2 spin model. Based on this as- 
sumption, we start from the mean-field Z2 planar state, 
only to derive the dynamical magnetic properties of this 
bond-type spin nematic phase. 

In the momentum representation, the Bogoliubov-de 
Gennes Hamiltonian for this Z2 planar state is given by 
^BdG=E.a4"^JfI°Vfe"and 



H 



(a)_JiD JiD 



(19) 



with Sf^ = sinfc^ and = cosfc^. The 4x4 7-matrices 
are defined as 



71 = (72 (X) (Tl = 



—iai 
iai 



(20) 



where the 2x2 Pauli matrices cr^ (/z = 1, 2, 3) in front of 
the (g)-mark is for the particle-hole space, while the other 
is for the spin space. Using the same notation, we define 
the other 4 anti-commutating 7-matrices as 72 = (12 cr2 , 

73 = 0'2 (Xi (T3, 74 = 0-3 (g) (To, 75 = CTi (g) Co- 

At this saddle point, the partition function is evaluated 

as 



Z^^\h] 
S^\h] 



exp 



1 ^ 

2-]^ETr(lnGoM/n), 



(21) 
(22) 
(23) 

(24) 



where the trace is over the momentum (fc), the Matsub- 
ara frequency (n), spin, and particle-hole indices. The 
single-particle Green function at the saddle point, Gq, is 
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given by 

^0,(fc.n,a|fc',«',a)W = ^n^n' Sk.k' Oo^k, lUJ^) 
1 1 ^ 

+ 2 /BIT' ^ ^ '5fc,fc'+q^ri,»i'+m /iq^^pM/i, (25) 
^ " /i=l q,m 

Qo^ik, iujn) = «w„7o - (26) 
where the bosonic fields are Fourier-transformed as 



laa 



}w=Jl[d-^'"'-'"'-'^t,i-) (27) 



with e„j = 2rmTf3 ^. The 4x4 Hermite matrices 
(/i = 1, 2, 3) are defined as follows: 

Ml = 7i5 = -ni75, «2 = 7i3 = -ni73, 

"3 = 735 = -«7375- (28) 



III. LOOP EXPANSION FOR CORRELATION 
FUNCTIONS 

In the previous section, we have introduced the Z2 pla- 
nar state as the saddle point solution of the general- 
quantum frustrated ferromagnetic model. Starting from 
this saddle point solution, we will calculate in this section 
the spin-spin correlation function, including the effect of 
the fluctuation of the auxiliary fields and the gauge fields 
around the local minimum within the random phase ap- 
proximation (RPA). Specifically, from the saddle point 
solution, we first obtain the Hartree-Fock contribution 
to the spin-spin correlation functions. We next give a 
symbolic expression for the spin-spin correlation func- 
tion, which includes the arbitrary order of the quantum 
corrections to this Hartree-Fock contribution. Based on 
this expression, we finally perform 1/iV expansion and 
calculate the first order quantum correction of 0{N), 
which corresponds to the random phase approximation 
term. 

Derivatives of the partition function Z[h] [Eq. ([13])] in 
external fields generate spin-spin correlation functionSfi^ 



with 



C^2(j,r) ^ {T{S^o':,msrAr)})\h=o 
N f 



(29) 



dK-^{Q)dh--^{T 



■ exp [ - NS] 



(30) 



/i=0 



Ar2 
dS 



dS 



exp [ - NS] 



(31) 
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FIG. 3: A contour-plot of (a) Imx£' (q, e) and (b) Imx+i (q, e) 
at J2 = l.lJi. The momentum q runs from (0,0) to (tt, tt), 
to (71,0) and back to the F-point (see Fig. (c)). Both of them 
have finite weight only at e > cJc, where ljc — 0.3. (d) Feyn- 
man diagram for ImxS]. 



where we do not take the summation over the flavor in- 
dex 'a'. The imaginary-time ordered correlation function, 
Cf^f^{j , r), is related to the real time dynamical suscepti- 
bility, Xw(9' = Xw(9' ^ e + iS), in the form 



E 

i 



r dTe-^^''-^+^--^C;;{j,T). (32) 
Jo 



If we omit the fluctuations around the saddle point 
solution, we obtain the Hartree-Fock contribution to the 
correlation function, 



-N- 



95(0) ggiO) 



h=0 



(33) 



h=0 
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Or equivalently, 



k,n 



Sq,odm,a 

16/3Na 



h.n 



Tr[go(fc',i'^«')'"i>. 



(34) 



k' .71' 



where the trace in Eq. ([M]) is over the spin index and 
particle-hole index. From Eqs. p6ll9l28p . it is clear that 
golkjiujn) is traceless for any /i — 1,2,3. Thus, we 
can omit the second term as 

TKT X! [doik + q, iujn + ie-m)Uf^ go{k, iuJn)u^] . 



8PNa 



(35) 



The Hartree-Fock solution Eq. (|35|) brings about only 
the so-called Stoner continuum, which stems from indi- 
vidual excitations of the 'free' neutral fermions (spinous). 
After the analytic continuation, ie„ e + iS, we obtain 
the real-time dynamical susceptibilities, whose imaginary 
parts are plotted as a function of q and e in Figs. [3] (a, b). 
The Stoner continuum thus obtained have a finite weight 
only above some critical energy, e > ujc- This feature is 
because the energy dispersions of the fermionic excita- 
tions, say ±Aj., are fully gapped in the Z2 planar state; 

A, = ^J',D'{sl + sD + Jlx^clcl + > 0-(36) 

One should also notice that continuum in Imxzziq, e) has 
no spectral weight at q = (tt, tt), because of the staggered 
[/(I) spin-rotational symmetry, Eq. (|T8l) . 

To capture the low-energy collective excitations, which 
emerge below these continuum spectra, we next include 
the effect of fluctuations of the auxiliary fields and the 
gauge fields around their saddle point values. The fluc- 
tuation fields for the Z2 planar state are given by the 
following 35 elements: 

rij,T) 

= [ReEx^3, lmEa,^3, ReEy^s, ImEy^s, 
ReD^^z, luiD^^s, ReDy,3, luiDy^s, 
ReEx,i, ImE^^i, ReEy^i, ImE'^^^i, 
ReDx^i — D, ImD^:,!, R^Dy i, ImDy 1, 
Rei?2;,2, Iini?:r,2, ReEy^2, Im£'j/,2 



ReDx 2, ImL':r,2, ReDy^2 - D, ImDy^2 
Rexx+y - X " 

1 ,-«2 ,„3 



- X, ^eXx-v - X, ^<2Vx+y - V, ^eT]x-y + V, 

Imxx+y, Imxx-y, Im772,+j,,Im77a;_y, ia]., ia^, ia^]. 

(37) 



Here, D^^^ and E^^^ with v = x,y stand for the /x-th 
component of the d- vector in the Cooper channel and the 
excitonic channel, respectively, both of which are defined 
on the nearest neighbor ferromagnetic links as 

whereas Xx±y/'>]x±y denotes the singlet 
Cooper/excitionic pairing defined on the x ± y-links as 

Vx±y(j + ^^^^^,t'^ = -i{fj,a[<J2]^pfj+e,±ey,p), 
Xx±y{j + ^^-^ 



In terms of these fluctuation flelds, the single-particle 
Green function appearing in Eq. ^ takes the form 

{k,n.a\k\n\a) 0,{k,n.a\k' ,n' ,a) 

{q,m)va{k,k'), 

V PjVa 



q m 



(38) 



where the bosonic fluctuation flelds r are transformed as 

The index a specifies the type of the fluctuation flelds 
(a = 1, • • • , 35). The summation over the repeated index 
a is made implicit and will be so henceforth. In the 
same sequence as in Eq. ([37]), the internal vertices Va are 
explicitly given by the 4x4 matrix forms 

v{k,k') 

= [-Cxl35, ~Sxll2, "Cy^a^, "Sy7l2, 
Sxfl: SxlU, Syfli Syfli, 
-CxflS, -Si, 723, -Cy7l5, -Sy723, 
-3x13, -Sa;734, -5^73, -Sy734, 
Cx731, Sa;725, Cy 731, Sy725, 
Sxl5, -Sxli5, Syj5, -Syfib, 

'-'a:+y74, ^a: — y74, ^2,^^72, ^x—yl'^, 
- Sx+y70, -sL-y70, C^+y724, C^-y724, "72, 724, 74] , 

(39) 



where 



Jl /^A" + 
— cos I — 

2 

Jl . fkfj^ + k' 
2 



J2 /' kx + k'^ ± ky ± k' 

Ox±y = YC0S( 

J2 . f kx ^ k'x i ky i ky 



s'x±y = ysm 



)■ 
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To obtain the quantum corrections to the Hartree-Fock 
contribution [Eqs. ([551) or ([55]) ]. we expand the action in 
term of the fluctuation fields r ; 



N 



a=l 

OO 

^ / , '-'ai.a2 



(40) 



71 — 3 



The coefficients in Eq. PO]) are given in terms of the 
internal vertex Va as 



c(") 

,a2 ,CKji ■ 



dro 



1 d^Si 



1 + 



2 dra^dra^ 



5n,: 



N 



with n > 1. The trace here is taken over the momentum, 
Matsubara frequency, spin, and particle-hole indices. 

To calculate the correlation functions Eqs. (15(11) and 
(|31]) . we expand the action in the small field with a 
certain flavor a in the form 

OO 

<,_^(0)^(2) ^^(n) 

^ 1 Oq, Q,/ / q / a' "T / ^ <-*ai,a2,--- ' cti ' (^2 ' 

n— 3 

-l 3 00 

//^l n— 1 
^ 3 00 

+ nT. E^2^;L-.".^"^---^->r^""' (42) 

/i.i/— 1 n— 

where the overhne stands for its evaluation at the zero 
field, e.g., = 5(o)U.o and t^^, = 5l'i,|,.=o. By 
definition, the expansion around the saddle point has no 

linear term in r. The coefficients 5^ ' ^ and 5^ ' ^ are 



given by 



.2 



/i=0 



2 dh^fdhl" 



h=0 



(43) 
(44) 



with ri > 0. In Eq. (jH]), we have also used the relation 

Z]fc,«Tr[5o(fc, «w„)i6^J = 0, which already ap- 
peared in the second term of Eq. The Hartree-Fock 

susceptibility xl^v'^^'' corresponds to S^i^ f^. 
The bilinear term in r is given by 



1 



1 



2 dvadra 



4:13N/ 



■'^Tr[gaVagaVa']. (45) 



Other 5^"^ (with n > 2), 5^^' 



, -?r(2,n) 

and o contain one 



loop composed by multiple one-particle green functions. 
For example, ' \ and ' ■* take the following forms 



(46) 



-^(2,1) 



16(/3A^a)3/ 



^ E {Tr[goM..goMA'9oVa] 



+ 



Tr[goWp9oMi/go^'a]}, (47) 



where k and iw„ denote the momentum and the fre- 
quency inside of the loop. The trace in Eqs. (|45l - 
I47p is only over the particle-hole and spin index. Ob- 
serving these expressions, notice that S^.'^^ ■■■ a„ and 

'^li^L^Qi,- - ,Q„ as well as s'^^^... ^^^ are generally functions 
of order unity in the large- TV limit. 

Based on Eq. spin-spin correlation functions 

= Xm",! + Xlt,ii can be expanded aai^ 



,.,l(g.»£n) = -- / 2^5^ „.„^...„^rai j 2^ ( 2^ '^ai,- ■ ■ ■ ^a, j 6 ° - ^ (48) 

n=0 m=0 ■ ;=3 

_^ „ OO OO 

1 1 1 



TO. 

m=0 j=3 



where iSq was omitted for simplicity. The Gaussian integrals over the real-valued fields r can be taken, reducing even 
numbers of the fields to a sum over all the possible pairwise contractions among the fields; 

J drri- ■■r2k exp[- NS^^^^,r^ra,] =J2 ^[^''^^^ \(i),<t(2) [5^^^' ^]^(3),^(4) ' ' ' [^^^^' \(2fe-i),<T(2fe) • (50) 
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FIG. 4: 1/A^-contributions to correlation functions, where a 
dotted line denotes external fields, a wavy line is the RPA 
propagator and a solid line is single-particle Green function, 
(a) contributes to while (b,c,d) to xji^,i(9, 

(a) spin-wave term, (b,d) Hartree-Fock (HF) term with renor- 
malized single-particle Green functions, (c) HF term with the 
vertex correction. 



The summation over a runs over the arbitrary permuta- 
tions among the 2k indices. 

To enumerate all the contractions possible in Eqs. (|48t - 
|49)) . use Feynman diagrams as depicted in Fig.|4l In these 

diagrams, all coefhcients S*^""* (with n > 3), S^^'^\ and 

S"^^'"^ in Eq. (|^ contain single loops made by solid lines 
and vertices. Each solid line denotes a single-particle 
Green function, and two neighboring solid lines are al- 
ways mediated by either an external vertex or an 
internal vertex v^. From the external vertex, a dotted 
line representing an external field is emitted. Any inter- 
nal vertex is connected with one and only one of others 
internal vertices by a wavy line representing the 'RPA 

propagator , [b \a^p- 

In the large-A^ limit, each single closed loop appear- 

ing in S (with n > 3) is in the order of unity due 
to the flavor degree of freedom and the prefactor 
Each loop connected to any external field (dotted line) 

appearing in S" or 5 ' is also in the order of unity 
as the flavor index of the internal line is set to 'a'. From 
Eq. (|5D)) . a RPA propagator is accompanied by factor 
N~^. This concludes that the leading contribution to 
the spin-spin correlation function in the large- iV limit is 
its saddle-point estimation, Eq. (PSI) [see Fig. [2Id)] and 
the first order quantum corrections in the order of 
are given by the Feynman diagrams depicted in Fig. 2] 
Those diagrams which vanish by themselves have been 
already omitted. 

Notice first that Figs. lUJb-d) take the same structure 
as that of the mean-field diagram depicted in Fig. Eld) , 
where the difference can be solely attributed to a proper 
renormalization of the single-particle Green function 
(b,d) or that of the external vertices (c). Thus, their ma- 
jor contribution is more or less modification of shape and 
intensity of the Stoner continuum. In Fig. Ufa) , on the 



other hand, the momentum and energy carried by one 
of the external lines are transmitted to the other only 
through the RPA propagator, into which various collec- 
tive excitations including magnetic Goldstone modes are 
encoded. Thus, some of low-energy poles of the RPA 
propagator generally show up as coherent bosonic peaks 
in the imaginary part of Fig. Ufa). 

To see this situation explicitly, we henceforth derive 
the expression for Fig. HJa). Notice first that a 4 x 4 
BdG Hamiltonian for the Z2 planar state satisfies the 
following symmetry relations 

73 iifc^ fc^73 - J^k^.ky+TT' 75 iifc^ fe^75 - fc^(51) 

These symmetries require that the RPA propagator de- 
fined in Eq. (j45p is always a block diagonal matrix with 
respect to the following four groups of fluctuation fields: 

ill =ReDy,3 el + ImD^.s e\ + Im£;^,3 Cg -I- ReS^.a e\, 

(52) 

R2 EEReD^,3 e\ + ImDj,.3 el + Im£'y,3 el + ReE,^^y, el, 

(53) 

Rz =ReL»j,,2 e\ + ReDy^i el + ImAj,2 63 
-I- lmD^,i el -f lTaEy^2 el + Im£':,.i 
+ ReE^^2 67 + ^eEy^i e| + Ivaxx+y eg 
-I- IvciXx^y e?o + ln\rix+y e\^ + Im77j,_y 2 
+ ial + ial 6^4, 

R4 EE(ReD^4 - D)e*+ {ReDy^2 - D) e\ 
^' ImAj.i 63 + ImD;r^2 64 + \vi\Ey,x e\ 
+ Im£;a;_2 Cg + Rci?2;_i 67 + Rei?y,2 e\ 
+ (R-eXa;+y - x) eg + (Rexa;-y - x) e^o 



(54) 



(Re772,+j^ - 77) -I- (Re?7a;_j, -I- 77) 642 

r^l3' (55) 



■ 2 4 

la^ e 



where {e^} denotes the orthonormal basis of 35- 
dimensional space R^^ and the renamed fluctuation flelds 
{Rp.,a} are defined through i?^ = R^^aS-^- Using this 
representation, the Gaussian part of the action is indeed 
decomposed into four parts, 



(56) 



where the matrix elements inside of each block are given 
as 



1 d^Si 



2 9r^,a9r^,^ 4/3 A^A 



-E 



X Tr[go(fe + q, iw„ + ?e„)'M^,/3go(fc,iw„)ti^,a] (57) 

in the momentum representation. The internal ver- 
tices here are also renamed such that "Y^^raVa = 
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One can also see that, under Eq. (|5T|) . any fluctua- 
tion fields in Eq. (|55p are disconnected from the external 
magnetic fields, while those in Eqs. (f52|) " (|54|) are coupled 
with the external vertices in the form 



(58) 



where 
^(1,1) 



4/3 A^A ^ 



r 

k' 



Tr [go(fe + q, iujn + if-n)Uf^ go{k, iw„)ti^ 



(59) 



Using Eqs. ((57)) and (|59)) . we finally obtain the contri- 
bution from Fig. IDJa) as 



^(14) 



^(1,1) 



Among this equation, S^.'^ clearly has no finite imagi- 
nary part, provided that |Rez| < minfc(Afc+q-|- A^), since 
Eq. ((59|) has the same structure as the Hartree-Fock so- 
lution Eq. p5p . Namely, for any fj, and a, it always has 
the form 



(Afc+q + Afc)2' 



(61) 



where the numerator is a regular function of z and real- 
valued on the real- axis of z. As such, any finite con- 
tribution to Imxjj^^j^ jj(<7, e) below the Stoner continuum, 
e < minfe(Afc+q -I- A^), can be solely attributed to the 
pole contributions in the RPA propagator, i.e., the zeros 
of the eigenvalues of Eq. ([57]) . 



IV. DYNAMICAL SPIN STRUCTURE 
FACTORS 

In this section, we discuss the dynamical spin structure 
obtained from the large- A'^ expansion up to the first order 

of A''"^. The spectra of ^"^nx^^zuil^ ^) ^'^'^ Imx+l 11(97 
at the zero temperature are calculated numerically based 
on Eqs. ([57 | - (|61l) . Their typical plots are shown in 
Figs. [S] and [5] for the parameter point J2/J1 = 1.1, which 
is inside the Z2 planar phase, i.e., Jd > 1.1 > Jc,2 (see 
Fig. [5]). For simplicity, the fluctuations of the temporal 
gauge fields, ia^ [ji = 1,2,3), are not included in these 
calculations, so that only the pairing fields comprise these 
collective excitations. 

In the next two subsections (Sees. ITV^ and HV Bl) . we 
first describe the character of these low energy collective 
modes, especially focusing on gapless modes at the mo- 
mentum points q = (0,0) and q = (tt, 0). When the 
relative ratio between FM Ji and AF J2 is changed to 
J2/J1 = Jc,2, the boundary to the neighboring U(l) pla- 
nar phase, two of the gapped bosonic modes at the (tt, tt) 



momentum point become gapless (see Fig. [7]), which cor- 
respond to the appearance of gapless gauge excitations 
('photon- like' excitation) in the U{1) planar state. In 
Sec. IIV Cl we will argue that, in the U{1) planar phase, 
the mass of two other gapped magnetic modes at the 
(tt, tt) momentum point simultaneously vanishes, leading 
to a breakdown of certain symmetries possessed by the 
Z2 planar state (see Fig. [S]). These symmetry breakings 
are actually consistent with the 'confinement effect' pos- 
tulated in this U{1) planar state. 



A. Near F-point 

The Z2 planar state breaks all the spin-rotational sym- 
metries, while preserves the translational symmetries of 
the square lattice. Thus, the RPA propagator generally 
has three spin wave modes which are gapless only at the 
F-point. The corresponding fiuctuation fields take the 
following forms at q = (0, 0): 



1 



V2 



(62) 



The 01-mode and 02,3-modes appear, respectively, as the 
gapless excitations in the longitudinal and transverse dy- 
namical spin structure factors (see Figs. [5] and [6]). Un- 
der the mirror reflection which exchanges the cc-axis and 
y-axis, the planar state is symmetric, and the (j)2- and 
(/>3-modes are interchanged. Thus, the 02- and 03-modes 
are energetically degenerate along (0,0) to {tt,tt), while 
the degeneracy is hfted along (0,0) to (7r,0) (see Fig. [S]). 

Near the F-point, the spectral weight of these spin wave 
modes always vanishes as a linear function of the momen- 
tum; 



e) = a\q\S{e - v\q\ 



(63) 



This is because the coupling between these spin wave 
modes and external fields are given by Eq. ([6T|) with the 
numerator of its right hand side vanishes as a linear func- 
tion of z and also because eigenvalues of the RPA propa- 
gator are always even in the Matsubara frequency. Under 
the Kramers-Kronig relation, the real part of the static 
susceptibility acquires a corresponding quantum correc- 
tion, which does not diverge at q = 0. This observation 
is consistent with the fact that broken spin rotational 
symmetry is induced by the ordering of the quadrupole 
moment. 

The vanishing spectral weight of the spin-wave modes 
at the F-point was also observed in other class of 
quadrupole order phasesji^ii^ and is a universal feature 
of spin-nematic states. A similar behavior of the dynam- 
ical spin structure factor is also expected in the collinear 
AF phase, which is stabilized in the strong AF J2 regime 
of the present J1-J2 model (J2 > 0.57Ji). To distinguish 
the current Z2 planar phase from this collinear AF phase, 
we also need to look into the spin structure factor near 
the (tt, 0) or (0,7r) momentum point (see Sec. IIVB|) . 
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FIG. 5: (Color online) Excitation energy spectrum in the dy- 
namical structure factor Imxzz(q,e) in the Z2 planar ground 
state for J2/J1 = l-l- The momentum runs along several 
high symmetric fc-points [see Fig. [3|c)]. The grey regions de- 
note Stoner continuum, (a) Spectral wight of the collective 
modes in Imxzz(q,e). (b) The momentum-energy dispersion 
relation, showing the collective modes given in the text, (c) 
The dispersion relation near the F-point is enlarged, showing 
the three gapped squashing modes. 



As the Z2 planar state also breaks the global SU{2) 
gauge symmetries, one might naively expect that the 
RPA propagator contains three other gapless modes, 
which correspond to the Bogoliubov sound modes. How- 
ever, when their couplings with 'gauge excitations' are 
included, all of these sound modes will be generally ab- 
sorbed into the gauge fields, only to endow the latter fields 
with finite Higgs mass:"^'^°'^^ The Z2 planar states do not 




(0,0) 




FIG. 6: (Color online) Excitation energy spectrum in the dy- 
namical structure factor lmx+-{q, e) in the Z2 planar ground 
state for J2/J1 — 1.1. The grey regions denote Stoner 
continuum, (a) Spectral wight of the collective modes in 

ImxH (<7, e). (b) The momentum-energy dispersion relations, 

showing collective modes and gapped squashing modes. 

have any physical Goldstone modes other than the spin- 
wave modes. 

The low-energy excitations below the Stoner contin- 
uum also contain well-defined massive modes near the 
F-point [see Figs. [5ljb,c) and Fig. Eljb)], which are direct 
analogue of the so-called squashing modes observed in 
the superfiuid "^He-B phasei^i^ They take the following 
form at q = (0, 0): 

pi = ai {el + el) + Pi {el + eg) + 71 (e?, + el^) , 
P2 = el- el, p3 = a2{el - e^) + ^^{el ~ ef), 
Pi = 64, P5 = 64, • • • . (64) 

Due to the SU{2) gauge symmetry, these modes are com- 
posites not only of the 'squashing' components in the 
Cooper channel, but also of those in the excitonic chan- 
nel. Note also that the /93-mode is a linear combination of 
the spin-wave component and the squashing component 
in the excitonic channel. This combination becomes pos- 
sible, because the spin-triplet excitonic pairing operator 
on the x-link and the spin-triplet Cooper pairing operator 
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on the y-link do not commute quantum-mechanically in 
the Z-i planar state (see Appendix). The spectral weight 
of all these gapped modes reduce to zero at the F-point, 
while some of them acquire finite weight for finite mo- 
mentum. In general, however, magnetic properties de- 
rived from dynamical spin structure factors cannot be 
a direct experimental probe for these squashing modes; 
their weight are quantitatively subdominant. 



B. Near (tt, 0)-point and (tt, 7r)-point 

When the momentum is away from the F-point, all the 
collective excitations are fully gapped, except for the ap- 
pearance of two gapless linear modes in Imxzz(<7,e) at 
q = (tTjO) and (0,7r) points [see Fig. [Sfb)]. The cor- 
responding fluctuation fields are given by Aq = eg and 
Cy, respectively. However, their spectral weight is essen- 
tially zero. This is because both fluctuation modes e| 
and 67 do not couple with the external magnetic field 

at all, \^t^\ci,z)\i = ^l^\q,z)\^ = at q ~ (7r,0) or 
(0, tt). This behavior should be contrasted to the dynam- 
ical spin structure expected in the neighboring coUinear 
antiferromagnetic phase, where the spectral weight re- 
mains finite even at the gapless points q — (tt, 0) and 
(0, tt), e.g. Imx;.;.((7r, 0) + 9, e) ~ h'8{e ~v!\q\) + -- - for 
\q\ < 1 and e < 1. 

Due to the mirror symmetry with respect to the (x+y)- 
axis, eigenmodes at the (tt, 7r)-point are always doubly 
degenerate. The lowest gapped excitations comprised of 
the fluctuation fields 



<P2 



bit 



7ei + 5(ei + e3o) 



(65) 
(66) 
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FIG. 7: (Color online) Dynamical structure factor Imxzz(q, e) 
at the zero temperature at the critical point between the 
Z2 planar and C/(l) planar phases, i.e., J2/J1 — Jc,2 = 
1.0448 (see text). The grey regions denote Stoner continuum. 

(a) Spectral wight of the collective modes in ImXzz(q,e). 

(b) The momentum-energy dispersion relations for the collec- 
tive modes, showing that the lowest excitations at the (tt, tt)- 
point become gapless. 



appear in ImXzz('Z, e)- Their momentum-energy disper- 
sion and the spectral weight near the (tt, 7r)-point are 
shown in Fig. [SJ The second lowest gapped modes com- 
prised of the fluctuations 

^z^eel + Cel, Lpi = ee\ + Ce\ (67) 

appear in Im^n (9:e)j whose energy dispersion and 

the spectral weight near the (tt, 7r)-point are shown in 
Fig. [6] When the momentum approaches this symmetric 
fc-point, the spectral weight of these four gapped modes 
vanishes as a quadratic function of the momentum, i.e. 
Imx^^((7r,7r)-hq,e) ~ a"\q\^5{e-m-v"\q\'^) for \q\ < 1. 
The vanishing of the spectral weight ai q — (tt, tt) is a 
consequence of the staggered C/(l) spin- rotational sym- 
metry in the Z2 planar state, Eq. (fT8| . As will be de- 
scribed in the next subsection, on decreasing J2, the ipi- 
mode and (p2-iiiode become gapless, only to comprise a 
photon-like dispersion at J2/J1 < Jc.2, where the other 
two gapped modes [1^93 and 954] simultaneously exhibit 
instabilities. 



C. Transition to the 17(1) planar state and its 
instability 

The saddle point solution suggests that, when the 2nd- 
neighbor antiferromagnetic exchange J2 decreases, the d- 
wave spin-singlet Cooper pairing amplitude 77 reduces to 
zero at the critical point J2/J1 = Jc.2 — 1-0448, while 
the spin-triplet Cooper pairing amplitude D and the s- 
wave excitonic pairing amplitude x remain finite even 
beyond the boundary, J2I J\ < Jc,2'^ Such a planar state 
in J2/J1 < Jc,2 becomes invariant under the following 
global [/(I) rotation around the 3-axis in the gauge space 

'^o^e-^^-'y^^''''^-^j (68) 

and hence that is dubbed as the U{1) planar state. 

Owing to the restoration of this continuous gauge sym- 
metry, some of gapped gauge bosons in the Z2 planar 
state become massless at J2I J\ = Jc.2- t/Ji-mode and (^2- 
modes defined in Eqs. (|65l66p are actually spatial com- 
ponents of these gauge bosons, which compose (part of) 
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a gauge-invariant Maxwell form for J < Jc.2 (see Fig.[7l). 
Specifically, a and 7 in Eqs. (|55|) and reduce to zero 
at J1I J\ — Jc,2i where these two modes become 



^e| + x(ei + e?o). 



Regarding that their amplitudes are infinitesimally small, 
one can identify these fluctuation fields as the spatial 
components of the gauge field introduced as 



T ^ 

3 



5]Tr[*].x^3e-'('^^-"?)-^vl,.+e.-eJ (69) 



with 



and 



The J7(l) local gauge 



symmetry implemented by Eq. (j68p requires that the 
spatial component thus introduced, in combination with 
the staggered component of the temporal gauge field, 
i.e. Oq = [—Xy^^iy ^, must take a gauge invariant 
quadratic form as their effective action, which turns out 
to be the Maxwell formiS 



gauge 



(70) 



The 'emergent' electromagnetic fields are defined as 



drtt^ — dattQ, and B = dxO^ 



dyal. 



In the 2 + 1 dimensional space, this Maxwell form does 
not suppress the fluctuations of these gauge fields effi- 
ciently, so that the U{1) planar state is generally un- 
stable against these fluctuations. That is, the space- 
time instanton which is allowed by the corresponding 
compact QED action, J (Pr{uE^ — K cos{eapdaa^p)}, 
proliferate in the 2-1-1 dimensional space, only to in- 
troduce strong confining potentials between two neutral 
'free' fermions (spinon)*22. In the context of spin-singlet 
quantum spin liquids, it is known that resulting confin- 
ing phases are accompanied by the reduction of the space 
group symmetry of original mean-field states.™ 

In the present situation, this symmetry reduction is 
driven by the condensation of the (pa-mode and ip^- 
mode defined in Eq. (j67l) . Namely, these two modes 
exhibit instabilities at the same time that the ipi- and 
(/?2-niodes become massless (see Fig. |5]). Since these 
bosons carry the (tt, 7r)-momentum, their instabilities 
break the translational symmetries of the square lat- 
tice (Fig. [IJc)). Moreover, their condensations also in- 
troduce finite z-components of the spin-triplet d-vectors 
connecting nearest neighbor sites (compare Eq. (j67p with 
Eqs. (|52l53p ). so that they also breaks the staggered U (1) 
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FIG. 8: (Color online) Dynamical structure factor 

ImxH (<7, e) at the zero temperature at the critical point 

J2/J1 = Jc.2 ~ 1.0448. The grey regions denote Stoner 
continuum, (a) Spectral wight of the collective modes in 

lmx-{ iQ,(-)- (b) The dispersion relations for the collective 

modes. At the (vr, 7r)-point, two collective modes become gap- 
less, only to condense at J2/J1 < Jc,2- (c) Staggered triplet- 
dimer state (left) and stripe- like order (right). 



spin- rotational symmetry defined in Eq. (IT51) . Having 
such 'triplet-dimer' orderings as the background, a pair 
of spinon and anti-spinon introduced on the uniform con- 
densate is expected to pay those energy cost which are 
proportional to the spatial distance between these two. 
Because of this strong confining potential, the pair are 
spatially confined to each other in the U{1) planar phase. 
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D. Transition to the tt-Hux states 

When the antiferromagnetic exchange J2 increases, the 
spin-triplet pairing field D gets smaller, while the other 
two remains almost constant. At J2/J1 > Jc,i — 1.325, 
the saddle point solution suggests that D = Q and 
X = 77 ^ 0, where the Stoner excitations become gap- 
less at five (inequivalent) symmetric momentum points, 
i.e. (0,0), (7r/2,7r/2), (tt,^), (tt, 0) and (0,7r). Corre- 
spondingly, all the collective excitations and their spec- 
tral weight in the spin structure factor merge into the 
lower edge of the Stoner continuum, when J2/J1 gets 
closer to the critical value Jc,i from below. 

V. SUMMARY AND DISCUSSION 

Using the standard large-N loop expansion, we have 
calculated the dynamical spin-structure factor of quan- 
tum spin nematic state in the spin-i square-lattice frus- 
trated ferromagnetic model. The saddle-point analysis 
in the large-N limit suggests that the 'd-wave' spin ne- 
matic state in this model is a certain kind of spin-triplet 
pairing states of the spinon field (neutral fermion), which 
can be referred to as the Z2 planar state. The spin struc- 
ture factors thus calculated have two aspects; they have 
both spin-liquid like character and 'spin-solid' like char- 
acter. The former feature is represented by the Stoner 
continuum which stems from the individual excitations 
of gapped free spinous. Owing to the ordering of the 
quadrupole moments, the dynamical spin structure fac- 
tors also acquire coherent peaks below this continuum, 
which signifies the existence of the gapless spin-wave 
modes. Being described as the a^m.-triplet pairing states, 
this state also has the so-called 'squashing' mode as its 
gapped collective excitations. All the gapless excitations 
including spin-wave modes have linear dispersions with 
respect to the momentum. Thus, temperature depen- 
dence of (the magnetic contributions of) the specific heat 
in the bond-type spin nematic ordered phase is given by 
a quadratic function of temperature, Cy ^ . 

Spectral weight of the spin-wave modes vanishes as 
a linear function of the momentum, when the momen- 
tum approaches the gapless F point. This behavior was 
also theoretically observed in the site-type spin nematic 
phase and presumably is a universal property of spin ne- 
matic states. Related to this behavior, Tsunetsugu and 
Arikawa previously calculated the temperature depen- 
dence of the longitudinal relaxation time of the nuclear 
magnetic resonance (NMR) in a site-type spin nematic 
phascjii To calculate the same quantity in the current 
bond- type spin nematic ordered phase, however, we need 
to estimate the 2-loop correction to the spin-spin cor- 
relation function in the present formulation. Namely, 
the most dominant spin-wave scattering process in the 
nuclear spin relaxation is usually the Raman process, in 
which a nuclear spin is fiipped by the simultaneous occur- 
rence of 1-magnon emission and 1-magnon absorption.^'* 




FIG. 9: Raman scattering process (2-loop level). 

Such a scattering process can be schematically shown in 
Fig. [9l which appears at the 2-loop level of the large- TV 
expansion. We leave this 2-loop calculation as a future 
open issue. 

The lowest gapped excitations at the (tt, 7r)-point are 
identified as a certain kind of Higgs bosons, whose fi- 
nite mass quantifies the stability of the Z2 planar state 
against the 'confinement effect'. Though being gauge-like 
excitations, these massive modes comprise not only the 
spin-singict pairing field but also the spin-triplet pair- 
ing fields. Owing to this triplet component, they have 
finite spectral weight even in the spin structure factor, 
once the momentum is deviated from the (tt, 7r)-point: 
the mass can be experimentally measured in terms of 
the inelastic neutron scattering experiment. When these 
Higgs bosons lose their mass, which is the case near the 
ferromagnetic phase boundary^i^ a 'linear' confining po- 
tential should be introduced between neutral fermions; 
the state is transformed into the other phase having no 
gapped free spinous. We found that, the moment this 
Higgs mass vanishes, the other two gapped bosons at the 
(tt, 7r)-point simultaneously change the sign of their mass 
from positive to negative, leading to their Bose-Einstein 
condensation. This condensation breaks the translational 
symmetry and a staggered C/(l) spin-rotational symme- 
try possessed by the Z2 planar state, which are consistent 
with the 'confinement effect' in the U{1) planar state. 
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Appendix A: on a gapped collective mode 

One of the gapped collective mode at the F-point is 
given by a linear combination of (gapless) spin-wave com- 
ponent and (gapped) squashing component; 

p^^a2{el~el)+l32{el-el). (Al) 

This mixing stems from the non-zero commutator be- 
tween the spin-triplet Cooper pairing operator and spin- 
triplet excitonic pairing operator. Namely, when defined 
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on two distinct links sharing one site, these operators do 
not commute with each other in general, 

Within the spin- wave approximation, the right hand side 
may be replaced by its expectation value with respect to a 
mean-field state. The Z2 planar state has a finite singlet 
Cooper pairing amplitude on the next nearest neighbor 
link. Accordingly, the diagonal components of the left 
hand side remain finite, whenever these two triplet pair- 
ing operators are defined on the nearest neighbor 2:-link 
and y-link sharing one site; 

[^{j-'-c,3),>y^E{3,]+y)A ^ i{fj-x,a['^-2\al3fj+y,l3) = ??(A2) 

Due to this commutation relation, the RPA propagator 
at the F-point has a finite coupling between the in-plane 
spin- wave component, ReZ?_ = RcZ)j;_2 — Rel?j,.i, and in- 
plane squashing component, Im£'_ = IvaE^ i — IvnEy 2] 



In fact, their off-diagonal matrix element ig is linear in 
the Matsubara frequency ie„ and linear in 77; 



/i 

/2 



8 iVA 



A. e^H-4Ar 



= /i 



4 TVa ^ 



si JlD'sl + Jlx'clcl 



4 TVa 



E 

k 



Owing to the point group symmetry of the mean-field 
state, these two components do not couple with other 
pairing fields at the F-point. Replacing ie„ by e + i5, one 
can readily see that one of the two eigenvalues is for the 
gapless spin-wave mode, while the other corresponds to 
the gapped collective mode; 



Im£'_ 



ReZ?- 



/i 









/2 



(A4) 







V 











(^3) 



The latter mode is nothing but the /93-mode defined in 
Eq. (jAip . where a2 is linear in 77. 
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